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In a generic gauge theory the gauge parameter dependence of individual Green functions 
is controlled by the Nielsen identities, which originate from an enlarged BRST symmetry. 
Qh! We give a practical introduction to the Nielsen identities of the Standard Model (SM) and 

to their renormalization and illustrate the power of this elegant formalism in the case of 
the problem of the definition of mass. We prove to all orders in perturbation theory the 
gauge-independence of the complex pole of the propagator for all physical fields of the 
SM, in the most general case with mixing and CP violation. At the amplitude level, the 
formalism provides an intuitive and general understanding of the gauge recombinations 
which makes it particularly useful at higher orders. We also include in an appendix the 
explicit expressions for the fermionic two-point functions in a generic gauge. 



1 Introduction 



Considering the subtle cancellations between various contributions necessary to make phys- 
ical observables gauge-parameter independent, it is not surprising that the variation of 
individual Green functions with respect to the gauge-fixing parameters are governed by 
symmetry relations. Formally, these relations can be shown to follow from an enlarged 
BRST symmetry in which the gauge parameters also undergo a BRST transformation (l|,|2[ . 
They are non- linear identities of the same kind of the Slavnov- Taylor Identities (STI), sat- 
isfied by Green functions at arbitrary external momenta, and are usually called Nielsen 
identities, after the seminal paper || in which they were first presented. 

The power of this technique lies in the possibility of factorizing the gauge parameter 
dependence in terms of new objects, the Green functions of BRST sources associated to 
the gauge parameters. In principle, this factorization holds at any order in perturbation 
theory, but its interplay with the renormalization procedure is not trivial. In the case 
of gauge-independent quantities, the gauge cancellations emerge from the recombination 
between these new objects and can be verified without an explicit evaluation of multi- 
loop diagrams. As we will see in the following, the mechanism of gauge recombination is 
revealed in great simplicity in the case of physical amplitudes. 

The Nielsen identities provide the appropriate framework to study any problem related 
to gauge dependence. They are well known to field theory experts and have been used in 
the study of the effective potential [D,|J and in high temperature field theory ||. Recently, 
they have also been studied in the context of the Abelian Higgs model |j| and of Yang- 
Mills theories with background fields. Our main purpose in this paper is to introduce 
the Nielsen identities of the full Standard Model (SM) and to study their renormalization. 
From a conceptual point of view, if all the physical parameters are fixed by normalization 
conditions directly based on physical observables, using the Nielsen identities it is possible 
to make sure that other quantities are gauge-independent. This applies to mass and mixing 
parameters of unstable fields — which we consider in detail — as well as to off-shell objects 
like effective charges. From a more practical point of view, we believe the Nielsen identities 
are also a useful tool for multi-loop calculations both in the electroweak SM and in QCD. 
The identities for the two, three, and four point functions that we obtain in the present 
paper can also be useful in this context. Throughout the paper, we will proceed in a 
pedagogical way and complement the formal treatment with explicit one-loop examples. 

As a demonstrative ground for the technique of the Nielsen identities we have chosen 
the problem of the definition of mass in the SM. This is an important and non-trivial 
issue which recently has received renewed attention pHllf, prompted in part by the high 



precision measurements of the Z mass at LEP and SLC It has been shown long ago |12 



that unstable particles are compatible with unitarity and causality. However, what makes 
the perturbative definition of the parameters associated to unstable fields a delicate and 
intriguing problem is the interplay between the phenomenon of resonance (which goes 
beyond perturbation theory as it implies the Dyson summation of an infinite number of 
diagrams) and the perturbative implementation of gauge symmetry. In particular, the 
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correct identification of the mass parameters of an unstable particle is not obvious: their 
gauge independence must be proved in full generality and their connection to experimental 
quantities clarified. 

A concept which is generally believed to be gauge-independent is the complex pole 
of the resummed propagator. To the best of our knowledge, there exists no general and 
rigorous proof that this is true. In this paper we use the Nielsen identities to provide the 
proof to all orders in perturbation theory and for all physical fields of the SM. The only 
assumption clearly needed to obtain this result is that the renormalization conditions for 
the physical parameters do not introduce spurious gauge dependence. This is the case 
whenever they are based on a well-defined set of physical observables. We also discuss 
how a mass parameter for the unstable fields can be consistently defined on this basis. 

We have organized the paper in the following way. In the next section we introduce the 
Nielsen Identity for the 1PI generating functional at the classical level. After a discussion 
of the renormalization conditions, we study the modifications introduced by quantum 
corrections in the most general scenario. In Sec. 3, as an elementary illustration, we 
discuss the Nielsen identities for the one-point Green functions. In Sec. 4 we consider 
the case of the W boson and prove the gauge-parameter independence of the complex 
pole of its propagator. Several comments and examples here should help clarify the most 
important points. As a digression, we also consider the infrared finiteness of the W pole 
mass. The analysis is then extended to the case of mixing. In Sec.[| we consider the 7, Z° 
sector and derive an interesting relation for the photon correlator at q 2 = in the SM. 
We then study in Sec. ^] the scalar sector and in Sec. |7] the fermionic sector. The following 
section is devoted to a discussion of the mechanism of gauge-cancellations in the case of 
four-fermion processes. Sec.|^ concludes the main text summarizing the most important 
points. We have collected some useful material in three appendices: in the first one we 
discuss some aspects of the derivation of the Nielsen identities and present the sector of 
the Lagrangian containing the BRST sources. In App. B we give a technical proof that is 
crucial for the results of Sec. 2. Finally, we provide in App. C the full one-loop fermionic 
self-energies in an arbitrary gauge. This completes the work of Ref. where the 
one-loop gauge dependence of the basic electroweak corrections has been considered. 



2 The Nielsen identities in the SM 

The idea behind the Nielsen identities is simple: the variation of the classical action with 
respect to a gauge parameter coincides with the BRST variation of a local polynomial 
in the fields. This is clearly necessary in order to guarantee the gauge-independence of 
physical observables. For example, the variation of an S-matrix element with respect to 
the gauge parameters corresponds to the insertion of the BRST variation of a local term 
between physical states, which is known to vanish. The Nielsen identities implement this 
simple idea at the quantum level. 

Our starting point is the Nielsen identity for the generating functional T at the classical 
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level @J 




(1) 



where \ = s£ is the BRST source associated to a generic gauge parameter £, s is the 
classical BRST generator, and Sr the linearized Slavnov- Taylor operator whose defini- 
tion is recalled in Appendix A. Notice that the operator coupled to x is non-linear in the 
quantum fields, therefore it requires a proper renormalization. The extended BRST au- 
tomatically takes into account the renormalization of the theory and the renormalization 
of the composite operators generated by the variations of the action with respect to the 
gauge parameters. Here and henceforth, we used the reduced functional, also defined in 
App. A, in place of the standard functional of proper functions. In the case of linear 
gauges, this allows us to write STI and Nielsen identities in a more compact way without 
modifying the Green functions of the physical fields. The 1PI Green functions of the the- 
ory are obtained differentiating T with respect to some of the SM fields. Differentiation 
of Eq. ([!]) therefore gives the gauge-dependence of a Green function in terms of products 
of other Green functions, which also contain the source x- 

We denote by r'"' (pi,j>2> •••) t ne 1PI Green function of cpi, cp 2 , ... at the ra-loop level. 
tpi can be any physical or unphysical field of the SM in a general covariant gauge, as 
well as any of the sources 7^, Xj associated to the BRST variation of tpi and of the 
gauge parameter T Vllp2 _ can be expressed as functional derivatives of the generating 
functional, the effective action T, with respect to the fields and sources •••> <p m , 



The exchange of two fermionic indices leads to a change in sign. We also adopt the short- 
hand notation for the partial derivative with respect to a generic gauge parameter £, 
whose associated source is generically called x- Some details concerning the action of the 
Slavnov- Taylor operator <Sr on T, the precise gauge-fixing term, and the complete source 
Lagrangian are given in App. A. Notice that V and its Green functions are renormalized 
objects, unless explicitly stated. 

Before we consider the quantum counterpart of Eq. ([]]), it is necessary to discuss the 
parametrization of the theory in some detail. We distinguish between three different 
categories of renormalization conditions. 

i) The ones that fix the physical parameters pi, namely the parameters of the classical 
gauge-invariant Lagrangian. They must be fixed using physical observables Oi (cross 
sections, decay rates, resonance parameters etc.): 



A set of renormalization conditions commonly used in precision calculations is given 
by the fine structure constant a, the Fermi constant Gp (measured in the muon 
decay), a s (measured e.g. from the ratio R of hadronic to leptonic decays of the Z°), 






(3) 
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the mixing parameters of the quark sector (measured e.g. in hadronic decays), and 
the masses of the Z°, the Higgs boson H, and all the fermions. In order to keep 
the renormalization program simple, it is indeed standard procedure to adopt mass 
parameters also for unstable fields. This has the advantage of establishing a direct 
connection between an experimental quantity and the two-point Green functions, 
T^q 2 ). On the other hand, the identification of the masses of unstable particles 
from the resonance parameters is not straightforward beyond the lowest orders of 
perturbation theory. For example, the masses of unstable particles are often defined 
in terms of the zero of the real part of the two-point function, i.e. by imposing 

Rer^(Mj) = 0. (4) 

This definition is not gauge independent beyond one-loop |8|,|10| unless is a stable 



field, but it has been used sometimes also in all-orders analyses [[14]] . A similar prob- 
lem of gauge-dependence may arise if one tries to define the mixing parameters in the 
quark sector in terms of two-point functions only, instead of relying on physical am- 



plitudes JT^]. It follows that, if physical amplitudes are not directly employed like in 
Eq. (^|), the consistency of the renormalization conditions has to be proved by means 
of the Nielsen identities, to all orders in perturbation theory, and the connection 
between theoretical constructs and experimental quantities has to be elucidated. In 
our discussion we will fix all the physical parameters using unambiguously defined 
physical observables (cross-sections, decays rates etc.). All sets of physical observ- 
ables are equivalent and are chosen according to the experimental precision of the 
inputs and to the convenience for the problem at hand. It is irrelevant for our analy- 
sis which set is actually employed. Having defined the physical parameters in terms 
of observables, we will show that the position of the complex pole of the propagator 
of all physical fields of the SM is a gauge independent quantity and can be used 
to define the mass parameters, provided the connection between the field-theory 
concept and the experimental quantities is clarified. This is the case, for instance, 
for the Z° mass parameter defined from the complex pole, whose relation with the 
resonance shape measured at LEP is well- understood — see the first of Refs. |§. 
This procedure applies also in the case of mixing between different fields. 

The conditions needed to prevent infrared (IR) divergences. Due to the presence of 
massless degrees of freedom, it is necessary to impose some auxiliary conditions that 
guarantee the correct IR behaviour of the theory. In particular, it is necessary to 
impose T AZ (0) = and similar conditions in the ghost sector [|14"1, |i~6|1 . 



iii) Other unphysical renormalization conditions, such as wave function renormaliza- 
tions, tadpole and gauge parameter renormalization. Apart from the case of the 
tadpole, discussed in Sec.|], we do not restrict ourselves to a specific choice, but 
simply require that they do not spoil the STI and do not affect the nilpotency of the 
Slavnov- Taylor operator. An alternative approach is followed in |14| . 



We recall that no invariant regularization is known for the SM. The implementation of 
dimensional regularization of Ref. M], for instance, is consistent but breaks the STI. They 



4 



have to be restored order by order through the introduction of non-invariant counterterms 



- see e.g. [[18|,[T9|]. This is a precondition to any discussion of the renormalization and it 
is necessary to recover the unitarity of the theory and the physical interpretation of the 
S-matrix amplitudes. 

Unlike the STI, Eq. ([I]) does not have to be preserved in the renormalization process, 
as the extended BRST symmetry is just a technical tool for the derivation of the Nielsen 
identities. Therefore, in the following we will consider the possible deformations of Eq. ([!]) 
induced by quantum effects in complete generality^] and write: 



dzT = Sr — T + A, (5) 




where the symmetry breaking term A is a dimension four operator with zero ghost number 
such that 5r A = 0. 

The investigation of the structure of A in Eq. (||) can be performed according to stan- 
dard cohomological techniques [T^,[T^,|TB|. Recalling that S£ = if S (T) = 0, the first 



step consists in writing A = X + iS r F with X ^ 5 r H. As can be intuitively understood, 
the part of A which can be expressed as the BRST variation of something else does not 
contribute to physical quantities. On the other hand, X does not decouple from the calcu- 
lation of observables and is usually called the cohomology of the operator «Sr- In the SM, 
X is composed of the dimension four gauge-invariant operators with zero ghost number, 
each of them representing a cohomology classQ The coefficients of the cohomology classes 
of Sr are the physical parameters of the theory. Therefore, a contribution to X can be 
absorbed into a renormalization of some of the physical parameters Pi and we can write 
X = Yli fii ■§-. r< For what concerns SrY, it admits different kinds of contributions and is 



extensively studied in the literature |lq, [21|,|22| . The most general expression for (|g) turns 
out to be 



s 5 r = (i + f) s T [jjL r) + ^-r + £ 7 « jv„ r + s,Jd 



In App. B we show how this structure is actually implemented and preserved at all orders. 
In Eq. (H) pi are the renormalized parameters of the SM, @f describes their explicit gauge 
dependence (or equivalently that of their corresponding counterterms), and ip is any of the 
physical or unphysical fields of the SM. When Eq. (^) is differentiated to obtain identities 
between Green functions, the operator J\f v counts the external fields, while p^, 7^ and 5 t 
parametrize the deformation of the Nielsen identity; they correspond to a renormalization 
of unphysical parameters. In particular, the third term in Eq. (^) renormalizes the external 
fields (wave function renormalization), the fourth renormalizes the tadpoles, and rescales 
the gauge parameters. As in the SM with restricted 't Hooft gauge-fixing there are four 



1 In the case of Yang-Mills theories, a discussion of the renormalization of the Nielsen Identity can be 
found in Ref.[2]; it agrees with the one given below. 

2 We recall that in the SM, besides the STI, some auxiliary constraints are needed to identify the gauge 
invariant operators. For a detailed discussion we refer to 
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gauge-fixing parameters & (i = Z, W, j,g) and as many sources Xh is i n f ac t a matrix. 
In the case of mixing between fields characterized by the same quantum numbers, 7* and 
Af v are also matrices. 

Eq. dH) shows the most general structure of the renormalized Nielsen identity. In many 
practical cases, however, the situation is simpler. First, our assumption on the renormal- 
ization of the physical parameters in terms of physical quantities implies automatically 
ft\ = 0. This follows directly from Eqs. (|3]), as Oi are gauge-independent physical objects. 
If the physical renormalization conditions were mistakenly chosen in a gauge- dependent 
way, non- vanishing (3f would arise because spurious gauge-dependence would be introduced 
inEq.d). 

In pure QCD, where naive dimensional regularization is consistent and respects the STI, 
it is customary to adopt a minimal subtraction (MS) as an intermediate renormalization 
condition^. Such MS subtraction leads in this case not only to /?f = 0, because the 
renormalized parameters are guaranteed to be gauge- independent P3| , but also to A = 0. 

One can also wonder whether the NI can be realized at all orders by an appropriate set 
of unphysical renormalization conditions. Although a complete analysis of this problem 
is beyond the scope of the present paper, the possibility of preserving the form of the NI 
without modifying ad hoc the structure of the ST operator (as in ||) seems unlikely ||. 
On the other hand, the point of view we have followed here has been to allow for arbitrary 
deformations of the NI. 

The decomposition of A in Eq. (|5]) into X and SrY becomes important in the calcu- 
lation of physical observables. Since any operator that can be expressed as the BRST 
variation of something else decouples from physical quantities, SrY is completely irrele- 
vant to their calculation. Hence, no contribution to the third, fourth and terms in the 
rhs of Eq. (||) has an effect on physical quantities. In Sec. |] we will consider, in particular, 
the gauge cancellations leading to gauge-independent physical amplitudes. Eq. (g) tells us 
that neither the renormalization of the fields, nor the one of the unphysical parameters, 
can spoil the gauge independence of the amplitudes. Only X — J2i Pi tt-T can make them 
gauge dependent [pq] . In other words, only the renormalization of the physical parameters 
of the theory affects the gauge-dependence of the physical observables. 



3 Tadpoles 

As a preliminary step in our analysis, we consider in this section the gauge-parameter 
dependence of the tadpoles. This is a very simple case and provides a first introduction 
to the use of the Nielsen identities. Because of the close connection between the mass and 
the tadpole renormalizations, the results of this section will be necessary in all subsequent 
applications. 

The 1PI generator T is defined as the Legendre transform of the connected generating 
functional Z. The Legendre transform is well-defined only if the linear terms in the fields 

3 This is also common in some one and two-loop electroweak calculations p3[. 
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(tadpoles) are removed at all orders in perturbation theory [PB|P^|. This is equivalent to 
setting the renormalization condition 

it^ = o, (7) 



and also corresponds to minimizing the effective potential at each order [28[ 



We now consider how the condition of Eq. (|7|) affects the Nielsen Identity. First, we 
differentiate both sides of Eq. (|6]) with respect to H. Taking into account Eq. ( |A3| ) and 
setting all deformation parameters to zero, but before employing Eq. (|7|), we obtain 

-$r K (o) = r X7HJ/ (o)r H (o) + r X7J/ (o)r HH (o). (8) 

All the external momenta are zero and we will drop them in the following of this section. 
As x is t ne source associated to a gauge parameter, it is a Grassman variable which does 
not depend on the space-time and does not carry any momentum. In deriving Eq. 
we have used the fact that the x's and the 7's have ghost number equal to +1 and — 1, 
respectively, and that non-vanishing Green functions must conserve the ghost charge. We 
have also used CP conservation to avoid, for instance, the appearance of H-G° mixing in 
higher orders. This assumption will be relaxed later. 

The renormalization of r x7H , which is logarithmically divergent, has to be fixed ex- 
plicitly. It follows from Eqs. (^) and (|7|) that for the Nielsen identity not to be deformed 
we must impose 

r<& = (9) 

at any order n of perturbation theory. If we allow the renormalization of T xlH to deform 
Eq. (H) according to Eq. @, however, we have 

(i + P «)r X7H + 5 i r HH (o) = o. (10) 

In the following we will consistently impose Eq. ([]). 

In the presence of CP violation, another tadpole amplitude emerges in the SM, con- 
nected to the vacuum expectation value of the CP-odd neutral would-be Goldstone boson, 
Gq. As the CP violation in the SM is confined to the fermionic sector, this will happen 
only at higher orders. In extended models, any neutral scalar field with zero ghost charge 
could develop a vacuum expectation value through radiative corrections. In all cases the 
linear terms in the fields must be removed. However, given Eq. (|7p, the STI imply the 
vanishing of tadpoles of the unphysical fields. Upon differentiation with respect to the 
neutral ghost field c z , Eq. ( |A3| ) yields 

= r c z 70 r Go + r c z^ H r H = o. (11) 

1,5=0 



SSr(T) 



Sc z (Q) 



4 Incidentally, it is interesting to see that the tadpole counterterm is generated by the BRST variation 
of a local counterterm: T CT — ST Sr (J d 4 xj H ) = ST J d 4 x gjpfe , where To is the tree level action and 
ST the coefficient of this counterterm. It then follows that a renormalization of the tadpole amplitude 
induces a shift proportional to ST in the mass parameters of all the SM fields. The previous equation 
uncovers also the unphysical nature of the renormalization of the tadpole. 
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Xw I i c^j 7ff 

** ~ ->-- ' 

Figure 1: One-loop diagram contributing to T xlH . 



To derive the previous equation, we have used Eqs. ( |A3| ) and (|A7[) and the fact that one- 
point functions are not vanishing only for neutral scalars with zero ghost number. As can 
be seen from Eq. (|A~7|), differs f rom zero already at the tree level, in which case it 
is proportional to v, the Higgs v.e.v.. From Eq. (|TTD it then follows that the vanishing of 



the CP-even tadpole implies the vanishing of the CP-odd tadpole at any order. 
Moreover, in the presence of CP violation a term S^ p J d 4 x 5T / '5G (x) should be added 
to Eq. (0). Using the STI for the two-point functions and the analogous of Eq. ([8]), and 
requiring 5f p = one then finds that Eq. (Qj is also valid, together with = 0. 

In the case of a model with two Higgs-doublets ||29||, Eq. ([TT]) takes the form 



r c 2 7o T Go + T c z lH T H + T c z lh T h + T c z 1A T A — 0. (12) 

where H, h and A are the physical neutral Higgs fields. It is sufficient to require only the 
vanishing of the tadpoles of the physical fields H, h, A. It then follows that the tadpole of 
the unphysical Goldstone boson is zero (identifying a flat direction in the Higgs potential 
PD|| ) at any order in perturbation theory. 

Before closing this section, it is instructive to check explicitly Eq. @ at the one-loop 
order. At this order none of the pathologies of (naive) dimensional regularization is ap- 
parent and we have an elementary example of a calculation with the x sources. To this 
end we expand Eq. (§) at 0(g) and consider dimensionally regularized Green functions 
before implementing the renormalization conditions. As a consequence of the Feynman 
rules given in App. A, the tree level Green functions T^ hH and vanish. We also have 

= by construction, while 1^1(0) = — M§. We therefore find 

d,r^=Mlr^ H , (is) 

where the last term is logarithmically divergent. It is straightforward to compute r*!| H 
using the Lagrangian given in App. A. Only diagrams of the kind displayed in Fig.Ol con- 
tribute and we recover the gauge dependence of T H given in Eqs. (11,12) of |13| . 



4 W boson 

As a first application of the technique to the case of the definition of mass, we consider 
the case of the charged W boson, which is particularly simple because it does not involve 
any mixing between different fields. We split the inverse W propagator into its transverse 
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and longitudinal parts 



r£W(g) = [gT - q ^f) KM) + q ^rKM)- (14) 

Our first aim is to obtain a Nielsen identity for the transverse part of the two-point 
function. The longitudinal part will be considered in Sec.|6]. As a first step, we differentiate 
both sides of Eq. ([!]) with respect to and W~ , take into account Eq. (|A3|) , and set to 
zero the Green functions which do not conserve the ghost charge. We obtain 

(q)T w ^(q))} (15) 



where t^ v = g^ v — q^q v jq 2 is the transverse projector and the superscript T indicates the 
transverse part of a Green function. From the discussion of the previous section we know 
that there is no non- vanishing one-point function and that T xlv for if = H,Gq, which 
describe the gauge- dependence of the tadpoles, must also vanish — see Eq. @ — if we 
impose 8% = 5f p = 0. The second line of Eq. fll5|), on the other hand, is not zero only for 
V 9 = Waj so that we obtain, at any order in perturbation theory (s = q 2 ^ 



d(Kw(s) = -2r T xlwW (s)T T ww (s), (16) 

w+w - = T l, w .w + - We now include tl 
using (3i = 0, Eq. (p!6|) becomes 



with = T^ +. We now include the possible deformations present in Eq. 



dF T ww {s) 



i + p^)r^( S ) + 7il KM- (17) 



For what concerns the mass parameter definition, the significance of Eq. (|T7|) is that a 
gauge invariant and self-consistent normalization condition on r^ w (s) can only be given 
at the location of the pole of the propagator. Defining the latter by 







we see that Eq. (|T7|) leads to d^T^ vw (s)\ s= s w = c^(r^ w (s w )), which in turn implies that 
the location s w of the complex pole of the propagator is gauge-independent at any order in 
perturbation theory. This is a remarkably non-trivial result of perturbation theory, as it 
concerns the parameters that describe the non-perturbative phenomenon of resonance. It 
relies exclusively on (3% — 0, which follows from our use of observables to fix all the physical 
parameters. The mass parameter m w and the width parameter Y w defined by s w = 
rr? w — i m w Y w are gauge independent quantities and, as a consequence of the discussion 
at point (i) in Sec. 2, m w can be adopted as renormalized W mass. Clearly, the precise 
connection between this parameter and related experimental quantities must be clarified 
in order to adopt Eq. fll8|) as a renormalization condition that directly fixes the W mass 
parameter. 
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Beyond one-loop order the definition of the mass parameter of an unstable particle 
in terms of its two-point function is not trivially gauge-independent as in Eq. (pi 



10]. Consider for instance the case in which the mass of the W boson is defined by a 



renormalization condition of the kind in Eq. (|4]), namely 

Rer^(A4) = 0; (19) 

the W mass counterterm is then ~ReT ww (M^). This is the conventional approach to one- 
loop mass renormalization |3TH33ll- Taking the real part of Eq. (|17]) at s = M^, expanding 
it at two-loop, and dropping and 7^ as they would not affect our conclusions being real, 
we obtain 

(9 t ReT T(2) = -2ReT T(1) ReT T(1) 4- 2Imr T(1) Imr T(1) 
t^rte j. ww z rte 1 mff at; l ww t a m± j. ww im ± xlwW 

where all terms are evaluated at q 2 = . Using the normalization condition Eq. ( |T9|) , we 



see that the last term is left over, so that Eq. (|PT|) is not satisfied. As a consequence, the 
mass parameter defined by Eq. ([0J) is gauge-parameter dependent beyond one- loop |L0| . 
As the imaginary part in the last term of the previous equation originates from gauge- 
dependent thresholds, there exists a class of gauges where Imr£^ w (M^) vanishes (Cf. 
Fig.^|) and for which the gauge parameter dependence of M w is only apparent at the 
three loop level ||10|| . The actual difference between the two mass definitions, AM 2 = 
Re[T ww (s w ) — T ww (Mw)], can be evaluated expanding T ww in powers of \s w — M^\ ~ 
T W M W = 0(g 2 ) up to 0(g 4 ). The result is AM 2 « M w r w lmT$w(M%r), which is clearly 
gauge parameter dependent. The renormalization condition (|l9|) is an example of definition 
of a physical parameter in a gauge-dependent way: beyond one-loop it induces ^ 0. 

A comment on the factor 7^ is now in order. As explained in the introduction, this 
term originates from the potential deformation of the Nielsen identity by the renormal- 
ization procedure. For instance, there is considerable freedom in the choice of both the 
wave function renormalization of the W field and the renormalization of T^, lwW (s). In 
case they do not respect the Nielsen identity, 7^ compensates for its breaking. Let us 
consider, for ex., the following two procedures at one-loop. A first possibility is to adopt 
a minimal subtraction (MS scheme) for both the wave function renormalization of the W 
and Tx-yvlw( s )- It should be clear that in this case 7^ = 0. A second possibility consists 
in using the on-shell scheme for the W field rescaling. If we now insist in using a minimal 
subtraction for T^ lwW (s), Eq. (|T6| ) is not satisfied by the finite parts of the counterterms, 
leading to a factor 7^ = ^x-ywwiM^)^ = |r^^(M^)|j^g, where the subscript MS 
means that only the finite part of this Green function is considered. Similar considerations 
apply to p£, which appears first at the two- loop level and is related to the renormalization 
of the gauge-fixing parameters. 

Like in the case of the tadpole, let us see explicitly what happens at the one loop 
level for regularized Green functions. Using Eq. flT3| ) and noting that the Green functions 
involving \ vanish at the tree level, Eq. ([15]) reduces to 



r^(s)+T^]=2TlfJ w (s)(s-Ml), (20) 
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Figure 2: One-loop diagrams contributing to Ty w j wW . 

where Tw is the contribution of the one-loop tadpole. The zero of the W inverse prop- 
agator is gauge-independent at s — M^,. Notice that T 1 ^^ w {s) describes the gauge- 
dependence of the residue of the physical pole, i.e. of the on-shell wave function renor- 
malization factor. An explicit calculation of the diagrams in FigfJ which contribute to 
^'xwlw w l ea d s t° the same ^-dependence of Aww reported in fl3fl ; the same happens for 
the £ Zi7 -dependence. 

We have seen that if the renormalization condition is not properly chosen, the mass 
parameter is gauge-dependent. A possible source of confusion, however, is the interplay 
of mass and tadpole renormalization. To make this point clear, it is sufficient to keep the 
discussion at the one-loop level. From Eq. fl2"0|) we know that the W mass counterterm 
5M^ = Rer^(M^,) + is gauge-independent. The tadpole renormalization according 
to Sec.|3], however, eliminates Tw from the previous expression and makes gauge- 
dependent. Nevertheless, we still have [3 M = 0. This is a consequence of the unphysical 
character of the tadpole renormalization. What is essential here is that the renormalization 
condition which fixes the physical parameter M w be gauge-independent, as is the case for 
Eq. ( |18D and not for Eq. ([19|). This and only this guarantees = 0. 

Two simple practical applications follow from Eq. flT7|), and we report them as illus- 
trations of the technique. First, we can consider the dependence of the W self-energy on 
the QCD gauge-parameter It is easy to show that the deformation parameters cannot 
affect it in this case, and that it is controlled by ^ -yww only. However, the ghost charge 
associated to the QCD gauge group and the one associated to the SU(2) group must be 
conserved independently of each other. Therefore, T^ slwW = at any order, which implies 
that the W two-point function does not depend on the gluon gauge parameter, as veri- 



fied in actual calculations at two and three loops ||34|| . The second application concerns 
the contributions to the W self-energy which are leading in an expansion in the heavy 
top quark mass. At the one- loop level, they are trivially gauge-independent, like all the 
fermionic contributions. At higher orders, one can use the fact that T^ wW (s) is only log- 
arithmically divergent to show that the gauge dependence of the heavy top expansion of 
T ww starts at the next-to-leading order. Again, this is not surprising, because the leading 
contributions in M t can be obtained in the framework of a Yukawa Lagrangian where the 
heavy fermions only couple to the Higgs boson and to the longitudinal components of the 



gauge bosons. This Lagrangian, which corresponds to the gaugeless limit of the SM |35 
does not require gauge-fixing. 
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Figure 3: Schwinger-Dyson equation for the W two-point function. The blobs on the inter- 
nal lines represent connected propagators (chains of bubbles), while the blob on the vertex 
represents a one-particle irreducible Green function. 

Infrared finiteness of the W mass 

The complex pole definition of mass based on Eq. ( |18D avoids also IR problems at higher 
orders in perturbation theory. It has been shown in Ref. |10[ that the use of the nor- 



malization condition of Eq. flTQ| ) leads to severe IR divergences in a class of higher order 
graphs containing the photon when the external momentum approaches the mass-shell of 
the W. As a consequence, in the resonance region, \s — M^\ <> M W T W , the perturbative 
series fails to converge, while it was found that the pole mass definition avoids all these 
pathologies. The origin of the problem is similar to the one of the gauge-dependence of 
the mass parameter defined by Eq. (|19D and is related to the need to take into account the 
imaginary part of T'^ vw in the renormalization procedure. 

More generally, the problem is common to all particles coupled to massless quanta, 
independently of whether they are stable or not, and concerns the perturbative description 
of the resonance region. For instance, in pure QCD it is well-known []36| that at two-loop 
order the two-point function of a massive quark is IR divergent at q 2 = m 2 unless the 
quark mass is renormalized on the pole. In Ref. JTTJ] it was shown that this property 
persists at all orders in QCD, namely that the perturbative pole mass of the quark in 
QCD is infrared safe (or finite). In the following we would like to approach the case of 
the W boson from a slightly different point of view, along the lines of ||11|| , generalizing 



some of the results of Ref. |T(J. We will show that the complex pole mass of the W is 
IR safe at any order in perturbation theory, namely that the renormalization condition of 
Eq. ([TJD does not lead to IR divergences in the resonance region of the W boson, nor to 
pathologies in the perturbative expansion. In that respect, the presence of the width does 
not alter the discussion in a relevant way. 

A convenient tool to analyze the IR properties of the W self-energy from a perturba- 
tive point of view are the renormalized Schwinger-Dyson equations (see e.g. [p6fl). These 
equations provide a simple iterative way to define the higher order graphs in terms of 
sub-diagrams. In the case of the W boson there are only two topologies containing the 
photon which should be considered, as they contain thresholds at s — and can lead at 
higher orders to IR problems. Their Schwinger-Dyson equations are graphically depicted 
in Fig.|3|. Diagrams with gauge-dependent threshold (like those with a charged Goldstone 
boson in place of the W) and with thresholds far away from the resonance region (like 
those with a Z° boson instead of the photon) can be discarded because their expansion 
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around s = M w does not contain non-analytic terms. 

We will treat explicitly only the case of the topology on the left side of Fig. [| as the 
other diagram can be discussed along the same lines. In this case the Schwinger-Dyson 
equation has the form 

r l 7 U- (p) ~ / dnk t %a pW , ^ p + k ) (*) z lu- ( k + p) r w +A 8W ; (*, P), (21) 

where r^ +w -(p) is the contribution to the self-energy due to the exchange of a single 
photon, T w + A w -(k,p + k) is the 1PI vertex, the superscript (0) indicates that the vertex 
is considered at the tree level, and finally Z A , g {k) and Z ( w+w -{k + p) are the connected 
propagator for the photon and for the W boson, respectively. To study the IR behavior 
of Eq. (|2l|) near the mass-shell, we now consider the transverse part of the self-energy 
T^\ w _(p) and approach the limit p 2 — > s w . We expand the propagator into the Dyson 
series of self-energies and tree propagators. Concerning the photon line, we recall that 
a convenient choice of the normalization conditions for the neutral gauge boson sector, 
i.e. r^ A (0) = 0, makes T^ A (0) vanish at all orders (Cf. next section). Therefore, the 
photon propagator Z AiyAr (k) is always proportional to 1/k 2 in the limit k — > and has IR 
dimension -2. 

For what concerns the W propagator, the IR divergent contributions are related only 
to the transverse component of Z w+W _(k + p) because the propagator of the longitudinal 
components of the W boson has a gauge dependent pole. In the on-shell limit for the 
momentum p and for k — > 0, the tree level W propagators present in the Dyson series 
for Z'^ r +W ^{k + p) are linearly divergent. Therefore, expanding Z W T +W ^ {k + p) around 
k = 0,p 2 = s w we have 

/ 1 \ n+1 

z c w l w -{k +p)\^- sw ~° E [ w - k ) lK w (Sw)] n • (22) 

Here we consider only the most dangerous terms, which vanish if and only if T' ww (s w ) = 0. 
Under this condition, Z^ + (k + p) is at most linearly divergent in the IR limit. 

w w p 2 =s w 

If, on the other hand, Eq. (^) is not satisfied, severe IR divergences appear in each order. 
The situation is not much improved if we move off the pole position in the resonance 
region. Indeed, in this case the W width acts as an IR regulator in the denominator of 
Eq. (p2|), but leads to a series where the denominator l/(s — s w ) ~ 0(1/ g 2 ) spoils the 
convergence of the perturbative expansion in the resonance region |10|| . 

The last information we need concerns the behavior of the vertex F' w+A ^ w .(k,p) (T 
refers to the transverse components of the W bosons) around k = 0,p 2 = M w . By 
analyticity and dimensional analysis, the vertex functions can be at most logarithmically 
divergent in the limit k — > (this can also be verified exploiting the STI together with a 
proper use of the renormalization conditions). Having IR dimension -3, it follows by power 
counting that Eq. (pT|) does not lead to IR divergences when the integral in the internal 
momentum k is performed around k = 0. 
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In summary, we have seen that the pole mass of the W boson, defined by Eq. (0), is 
IR safe to all orders in perturbation theory and that only if this definition is adopted a 
perturbative description of the resonance region is possible. 



5 The Z — 7 system 

The main difference between the case of the W boson and the one of the neutral vector 
bosons is the presence of mixing. We now directly use Eq. (|>D with /?| = and set = 
for ease of notation (doing otherwise would not modify our results). Following the same 
steps as in the derivation of Eq. fll7D , and keeping in mind that the abelian vector field 
does not need a BRST source, we find for i,j = A, Z 



c w TUs) ~ s w Tf A (s)} (r T xr Js) - 7 y , (23) 

where t| 3 = 7^ is the deformation induced by the possible mismatch between the wave 
function renormalization matrix Zij and the renormalization of T xi3 j. We recall that 
rj k (s) is a symmetric matrix. We now consider the quantity 

V T AZ {*) = det ( ll A i S } If (24) 



which appears in the denominator of the propagators of the photon-Z system (see for 
ex. |33[1). If we are interested in the analytic structure of neutral current amplitudes in 
the typical configuration of a high-energy collider, where external fermion masses can be 
neglected, V AZ (s) is what we need to investigate. It is straightforward to derive 

d{D T AZ {s) = -2 (c w T T xyiZ {s) - c wl \ z - s w T T xj3A (s) - s wl l) V T AZ {s). (25) 

This tells us that the zeros of V AZ identify gauge-independent quantities. On the other 
hand, we know from the STI that V AZ (0) = (see for ex. |32| ; Ref. [IB| considers also 



the case of CP violation) which in turn implies by analyticity V AZ (0) = 0. This result 
ensures the existence of a massless state, the photon. T> AZ has, however, another zero, 
corresponding to the Z° complex pole, at q 2 = s z . As in the case of the W boson, this 
result implies that the position of the complex pole is a gauge independent quantity and 
that the only self-consistent normalization condition for the Z° mass is the one given in 
analogy to Eq. flTg|). With the exception of the IR problems, all the discussion on the W 
mass applies directly to the case of the Z° boson mass || . A Ward Identity similar to the 
Nielsen identity of Eq. ( |17D has been applied in to the case of the Z° resonance, to the 
same avail. 

Another interesting application of Eq. fl23[) concerns the photon correlator at q 2 = 0. As 



is well known [ 3^| , using the renormalization condition T^ z (0) = the result T> AZ (0) that 



we have used above implies T^ A (0) = 0. In this case it is straightforward to verify from 
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Eq. ( P3D that the derivative wrt q 2 of the photon two-point function calculated at q 2 = is 
gauge-independent at all orders. Imposing the condition T^ z (0) = in the expression of 
c^r^ z (0), we obtain the constraint r^ 73A (0) — 7I4 = 0. We can now differentiate d^T T AA wrt 
s and evaluate it at s = 0. Using the various constraints we have obtained, we immediately 
derive 

d 



= 0. (26) 

s=0 



Notice that no particular renormalization condition on the derivative J^r^J s= o has been 
imposed, so one should think, for instance, of a minimal subtraction. This interesting and 
non-trivial result shows that under the condition T^ z (0) = and at s = there exists in the 
full SM something analogous to what happens in QED, where the vacuum polarization of 
the photon is gauge- independent for any s (see for ex. Ref. ||37|| ). An alternative derivation 
of Eq. (pH) can be obtained starting from the physical photon-electron amplitude at s — 0, 
proceeding along the lines of the discussion of Sec.|§ and taking the gauge-independence 
of the on-shell amplitude for granted. 



6 The scalar sector 

In the previous section we have studied a first example of mixing. Indeed, mixing occurs 
in several other cases in the SM and in most of its extensions; all can be treated in a 
way very similar to the {7, Z} case discussed above. In this section, we first consider the 
matrix T^(s) of the two point functions relative to the scalar fields <fi = {0i, 02, 0n} i n 
the general case of mixing and show that the gauge dependence of its determinant follows 
an equation analogous to Eq. fl2"5|) , if the rank of T^(s) is equal to its dimension n. As CP 
violation is present in the SM, we then consider the system formed by {A L , Z L ,G , H}, 
where the subscript L denotes the longitudinal component of the vector boson fields. This 
system is highly constrained by the STI and we show that in this case the complex pole 
of the only physical field, the Higgs boson, is gauge- invariant. In an analogous way one 
can consider the {W^, G^} system, which however has no physical degree of freedom and 
is completely constrained by the STI. 

The general form of the Nielsen identity in the case of a system <p of fields characterized 
by the same conserved quantum numbers can be obtained in analogy to Eq. (p3|) and reads 

d ( T*(s) = A(s) T*(s) + T+(s) A'(s), (27) 

where we do not need to specify the matrices A and A' any further. Using In det = 
tr In and exploiting the properties of the trace, one finds for T>§ = det T* 

d(D^s) = tr [A(s) + A»] V^s), (28) 

which generalizes Eq. (]25f) in the case the rank of T^(s) at arbitrary s is equal to its dimen- 
sionality. In the case of n scalar fields this ensures the gauge-independence of n complex 
poles. Notice that the physical information contained in the matrix is not restricted 
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to the physical poles. Indeed, the higher order definition of the mixing parameters is af- 
fected by the off-diagonal elements of T^. In general, it does not seem possible to form 
gauge-independent quantities on the basis of two-point functions only, i.e. of T^, and to 
employ them to renormalize the mixing parameters On the other hand, the mixing 
parameters can be safely defined in terms of physical observables such as mesonic decay 
rates. 

Neutral current processes are mediated by photons and Z°, as well as by scalar fields, 
like Go and the physical Higgs. As it is well-known, the propagator matrix is obtained by 
inversion of the two-point function matrix and, in the process of inversion, the transverse 
and longitudinal components of the vector boson fields decouple. Having considered the 
transverse degrees of freedom in the preceding section, we can now limit ourselves to the 
system formed by the longitudinal components of the photon and of the Z° and by the 
Higgs and the neutral Goldstone bosons, which we denote by S = {A L , Z l ,Gq, H}. The 
two point functions involving one vector boson and one scalar are defined extracting g M . 
In this way, T is the 4x4 matrix of the two-point functions of S. 

The system S includes unphysical degrees of freedom. As we have noted in the intro- 
duction, even at the tree level the Green functions of unphysical fields are modified by 
the choice to use the reduced generating functional T in place of the complete functional 
r c (see the App.A). For the purposes of this section, however, the reduced functional 
simplifies significantly the derivation without affecting the physical information we can 
extract from T . In a way, this can be viewed as a consequence of the fact that the cancel- 
lation between the unphysical degrees of freedom occurs independently of the gauge-fixing 



sector 



Each row of T s is connected by a STL For instance, differentiating Eq. ( |A3| ) with 
respect to and c A , we obtain for the first row 

c w s w^cais^) r aa i s w c w "I - c w V CA ^ 13 ) T AZ + T AGo T CA ^ + Y AH T CA ~ fH = 0. (29) 

Similar identities can be derived for the other rows, so that the STI for the two-point func- 
tions can be written as T S V CA = 0, where V? A = (c 2 w -s w T CAli , s w c w +c w T CAl3 , T CA7o , T CAlH ). 
Since <fi includes the unphysical components of the photon and Z° fields and since we have 
eliminated the gauge fixing sector of the tree level Lagrangian in using the reduced func- 
tional (see Eq. (|A2|) ), it is perhaps not surprising that there is no propagator for A L and 
Z L and that detr s = or the rank of T s is less than 4. In fact, T s has another lin- 
early independent eigenvector V Cz with zero eigenvalue, corresponding to the set of STI 
obtained by differentiation wrt c z . Therefore, the rank of T s is at most 2 and that we 
cannot use Eq. (f2^) at this stage. Moreover, the sub-matrix of F s identified by the indices 
Go and H has the same rank as T s . This can be seen by noting that V Cz and V CA can be 
orthogonalized in the subspace of the A L and Z L components because 

det ( ^ ~ SwT ^_fz s ^ ~ ^r cz73 \ =1 + o ^^ 0} (30) 

\ S W C W + C W L c A -y 3 S w + C W L c z ~i3 ) 

Having eliminated the unphysical longitudinal components of the vector bosons, we can 
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now concentrate on the sub-matrix 



YH _ 1 G Go 
V ^HG 




(31) 



whose rank is equal to the one of V s . Indeed, at arbitrary q 2 , its rank is 2, so that 
Eq. ( p8[) is satisfied. T n is very similar to the 7 — Z° transverse mixing matrix. Even 
if the CP violation mixes up physical and unphysical scalar fields at high perturbative 
orders, it is not difficult to disentangle them taking advantage of the STL At q 2 = the 
two STI obtained by differentiating wrt c A)Z and G imply that detT w (0) = 0. This zero 
is related to the G° field and is located at q 2 = (in the standard gauge it would be at 
q 2 = i z M 2 ) as a consequence of the use of the reduced functional. The remaining zero, at 
q 2 = s H , corresponds instead to the physical pole of the Higgs boson and its location in the 
complex plane is therefore gauge- independent, as it follows from Eq. fl28|). A discussion of 
the relation between the pole mass and the conventionally renormalized mass of the Higgs 
boson in this case can be found in Ref. 



7 Fermions 

The treatment of the fermionic sector is only slightly more involved than that of the scalar 
sector. Again, we consider the most general case of mixing and call l?-' the matrix of the 
fermionic two-point functions, IVf- In the case of massless neutrinos, there is no mixing 
in the leptonic sector and T lept is a diagonal matrix. As a first step, we need to decompose 
into scalar pieces: 

r'(p) = S L (p 2 ) p!P L + E R (p 2 ) i>P R + Z D (p 2 )P L + ZUp 2 )Pr (32) 

where P LyR = |(1 =F 7s) are the left and right-handed projectors. As can be seen by 
inverting , the relevant quantity for the fermion propagator matrix is the matrix |[{J 

K f (p 2 )=p 2 i; L -^ D Z- 1 Z D , (33) 

where we have dropped the p 2 dependence of the £ matrices. Since the determinant of this 
matrix appears in the denominator of the fermion propagators, we want to study its zeros, 
i.e. the zeros of the eigenvalues of Kj. We recall that by pseudo-hermiticity T$ = 7q 7 , 
so that T^ L {p 2 ) = T, L (p 2 ) and T> R (p 2 ) = T* R (p 2 ) (this is actually true below thresholds, but 
it does not affect our conclusions). Hence, the matrix Kf(p 2 ) is hermitian and can be 
diagonalized by means of a unitary transformation. Under the usual assumption flf — 0, 
the gauge-parameter dependence of T^ is described by a Nielsen identity which has exactly 
the same form of Eq. (|2"T|). Setting furthermore = 7^ = for ease of notation (the results 
would not change), we have A = —T x j v ^ and A' = — T xf j f p, which have a Dirac structure 
and undergo a decomposition analogous to Eq. (^). Again by pseudo-hermiticity, we find 
that in this case A L fl = (A' L R Y and A D = A' D . It is then straightforward to verify that the 
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components of satisfy 



d^ R = A D Z R + E fl At + A R tf D + Z D A R (34) 
d^ D = p 2 (A R Z L + Z R A[) + A D Y, D + Y, D A L 



from which it follows that 

d ( K f = K f F + F ] K f , (35) 

with F = A D — A^ 1 ^. Without using pseudo-hermiticity, we would have F' 7^ F* 
in place of F^ in the previous equation. Eq. fl35j) is in the form of Eq. (|27f) and therefore 
T>f = det Kf satisfies Eq. ( p£f) . We have therefore algebraically reduced the problem in 
the fermionic case to the scalar one. In the case of mixing between n fermions, the gauge- 
parameter independence of n complex poles is thus warranted. Again, this result holds for 
any choice of the fermion wave function renormalization and relies solely on the flf = 
assumption. 

The above proof is new and valid in the full SM. For what concerns pure QED and QCD, 
the result that the pole masses of the electron and of the quark are gauge-independent is 
not new and has been obtained both using the Nielsen identities |37] and in different ways 



TT. JD|. In QED (QCD) the situation simplifies considerably: writing Tjj = B $ — mA, 
where m is the mass of the electron (quark), and decomposing A in an analogous way, we 
find 

2 

— dcA = — B A B + mAA A ; -d*B = m {A A B + B A A ) , (36) 
m 



which could be tested up to 0{a 2 s ) against the general gauge calculation of Ref. 



The proof of the IR finiteness of the fermions in the SM follows Ref. jjlTJ and the final 
discussion in Sec. 3 and is already present in nuce in Ref. [0]. For completeness, in App. B 



we present the explicit gauge-parameter dependence of the one-loop fermionic self-energies 
in a general gauge for the full SM. Remembering that A first occur at the one-loop level, 
it is straightforward to see that they satisfy the Nielsen identities Eq. (|34]). This completes 
the set of expressions given in Ref. [[HJ and is very useful in particular applications. For 
instance, Eqs. ( Pl"l - |C3| ) have been used in Ref. [[DJ to discuss the gauge dependence of 
the one-loop definition of the CKM matrix. Indeed, as noted in the previous section, 
the renormalization of the mixing parameters is a delicate subject for what concerns the 
gauge-parameter dependence. An adequate framework for studying it is the Background 
Field Method Hl6| . In the case of the fermion mixing a comprehensive analysis has been 
presented in Ref. \T5\. 



8 Application to physical amplitudes 

In this section we apply the formalism of the Nielsen identities to four-fermion physical 
amplitudes and study the mechanism of gauge cancellations at any order in perturbation 
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theory. Our purpose here is not to prove the gauge-independence of the physical am- 
plitudes, a result which was accomplished in full generality long ago at the level of the 
generating functional ||25|| . We would rather like to study a specific example and carry 
out the analysis at an arbitrary order in perturbation theory. The use of the Nielsen 
identities allows us to uncover the regularities of the gauge recombinations between the 
different components (vertices, boxes and self-energies) in great generality. The following 
derivation is formally independent of the perturbative expansion of the Green functions. 
In other words, if we work at order n in perturbation theory the Green functions have 
to be expanded up to this order, but the factorization works independently of that. At 
the one-loop level, a similar factorization is also accomplished diagrammatically by the 
Pinch Technique (PT) j|2| , whose extension at higher orders has however proved problem- 
atic. Unlike the PT, the Nielsen identities control only the gauge parameter variation and 
cannot be used to construct explicitly gauge-independent proper functions which satisfy 
basic requirements and tree-level-like Ward identities. However, they may prove useful 
in the search for the higher-order extension of the PT. The analysis of this section gives 
us also the opportunity to present explicitly the Nielsen identities for vertices and boxes 
involving fermions, which are interesting in their own respect as they appear in most 
phenomenological applications. 

We first consider the truncated Green function Zfj^ c N (see e.g. JZf]]) for a generic four 
fermion process fjfj — ► fgfN an d we decompose it in terms of irreducible diagrams and 
propagators. We will use capital and lowercase letters to denote fermions and bosonic 
fields (scalar as well as gauge vector bosons), respectively. Therefore, Zfj and ZL are the 
propagator functions of fermions and bosons. Following the convention of the preceding 
sections, irreducible boxes and vertices are denoted by Tjjkni ^Iju aric ^ ^jkn- To keep the 
notation simple, we drop Lorentz indices and the dependence on the external momenta. 
The physical amplitude M.jjkn f° r our process is obtained from Zf^ c N using the LSZ 
reduction formula ||26|| , which in the case of fermion with mixing reads |32[ 



M IJKN — lim Zjj, Zj,j Zpj?£, N , Z Rk ,Z^, n , (37) 

on— shell 



where the on-shell limit includes the projection on the asymptotic states and Z controls 
the relation between the asymptotic states and the renormalized spinors: 

Z) l jU as j = u I . (38) 

The matrix Z can be computed from the conditions |32| (quantum equations of motion) 



Tjj uj(mj) = 0; w/(m/) Tjj = 

Zjj Z 

" TjjUj{mj) = uj(mi); u I (m I )T n =«i(mj), (39) 



rrij p — mi 



using the fact that (tf — rrij) u as {m j) = at any order by definition. Of course, Z should 
be decomposed in left and right-handed parts, Z = Z L P L + Z R P R . Notice that the first line 
of Eqs. (^9]) implies det Kf = and consequently includes the requirement that the mass 
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4>'j # 




Figure 4: Nielsen identity for the two-point function Y^.. 

parameters of the external fermions are renormalized on the poles of the propagators (see 
Sec.0). Strictly speaking, the LSZ formalism applies only to stable external states, i.e. to 
the electron and neutrinos and, to a good approximation, to the muon. Nevertheless, we 
will consider here the general case of mixing. We also stress that the LSZ factors Z should 
not be confused with the wave-function renormalization factors for the external fields. 
Of course, the latter can be chosen by imposing Eqs. fl39|) together with Z = 1 (on-shell 
scheme ||32|| ) , but there is in general no restriction on their choice (see also [pl|]) and it 



is even possible to avoid them altogether, in which case Z is divergent. Once the wave- 
function renormalization has been defined, for instance through a minimal subtraction, 
the factors Z can be computed from Eqs. 



As a first step, we consider the gauge variation of the truncated Green function Zj™^ 
In the most general case of mixing, Z l f^ c N is decomposed in the following blocks (we sum 
over repeated indices) 

ZfjRN — i^ijRN ~ T/jj Z^ TjK N — Tj Ni Z^j r^j, (40) 
from which we obtain 

- {d^iJi) %ij ^jRN - T/jj (ptZfj) ^jkn - T/jj Z\j d^Tj^ N (41) 

We can compute the different contributions d^Y jj N j^,d^Y and d^Zfj using the Nielsen 
identities. The identity for the propagator functions Zf- and Zjj is easily derived from 
the identity for the irreducible two-point functions Yij and Yjj. As we have seen, the 
general form of the latter is 

ab Yy^i c r^, a ^y c Tcb, (42) 

where the indices a, h apply to both the bosonic and fermionic case. As usual, we employ 
the procedure of Sec. 2, remove all tadpoles, set 5 t = <5f p = 0, and assume 4 = 0, Vz (this 
is consistent with the LSZ use of pole masses, as we have seen). Concerning the and 7^ 
factors, we avoid them here in order to keep the formulas simple. However, following the 
discussion in Sec. 2, they are bound to drop out of the amplitude and this can be explicitly 
verified. Eq. (|42]) can be graphically represented in the very simple way shown in Fig.^j. 
Notice that the momentum flows along the horizontal line and that the insertion of the 
static source \ does not carry momentum, unlike the one of r ) tp . 
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Using the relations Z?- Tj k = i8ik and Zj K T^j = iSijl, where 1 is the identity matrix 
for the Dirac indices, we obtain the Nielsen identities for the propagator functions, 
which read 

fy Z ij = Z ik^xkjj +^xnk z ip (43) 
d £ z h = z lK T x i<vj + T xviK z< kj, (44) 

for bosons and fermions, respectively. Graphically, these identities can be represented by 
Fig. [| after replacing the blobs with the \ insertion by their mirror images and exchanging 
the corresponding indices. For the three-point functions we have 

+ ^IKi^xmJ + ^IK^xiVKJ + ^xIVK^- KJi + ^xilVK^KJ • (^5) 

We see that the gauge-dependent terms of the form of F xyi j introduced by the propagators 
in Eq. ( f4"T| ) are exactly cancelled by the last term in the first line of Eq. (|4"5|), i.e. by the 
vertices alone. Therefore, the boxes are not necessary to remove the gauge-dependence of 
the internal self-energies. The identity for the four-point functions is 

~ 9^ijkn — r/ Jm r x7m ^jv + Tj Nm T xlm Kj + T xlm jj T m x N + T xlm j N r m Kj 

+ r xf?s jk N + r x ijfj S N + ^ x jjk vs Fg N + t x j vs ^ n t§j (46) 

+ ^I.JKS ^XVsN + ^ISKN ^XVSJ + ^xhs ^SJKN + ^x^VS ^IJSN ■ 

We now distinguish between the different Green functions containing the source \: 

1. Terms of the form T xl jj are present both in the gauge variation of the boxes (first 
line) and in the one of the vertices (first term). They cancel against each other in 
the sum ( ffi"D according to the pattern 

r x7m jj r, m j z\- Tjk N — i r x7m /j T m ^ N = o , 

-{Win) % TjKN ids^UKN 

where we have specified which part of Eq. fl4l] ) generates each term. 

2. The factors containing T xifjK j in the second line of Eq. ( f45|) and T xfjsJ g N (the whole 
second line of Eq. ([IB"!) ) always multiply a two-point function of the external fermions 
like Tjj. When they are contracted with the external spinors, these terms vanish, 
as a consequence of Eq. . 

3. The remaining terms contain Green functions of the kind T X f II j and J^ x j^j which 
multiply vertices and boxes in Eqs. (f45|) and fl46|) , respectively. As we will see in a 
moment, they are cancelled by the LSZ factors. 

Adding together the various pieces, the gauge-parameter variation of the on-shell truncated 
Green function can be expressed in terms of the truncated function itself: 

a ytrunc I t-i ytrunc _i_ "p _ ytrunc , ytrunc t-i , ytrunc p (a n\ 

~ u £, Zj IJKN\on-shell — *■ x In s ^SJKN <~ L x^vs "IJSN ' ^ISKN L XVsJ *+" ^IJKS 1 XVsN, < ) 
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according to the usual form for the Nielsen identities. Of course, this on-shell factor- 
ization holds in general for any amputated Green function, as it follows from the gauge 
independence of the S-matrix. 

We are now ready to apply the LSZ reduction formula. The gauge variation of the 
factor Z can be computed from Eq. ([38]) and Eq. (|39D using the Nielsen identities for the 
two-point functions and the gauge- independence of the asymptotic spinors u as j. We then 
obtain 

lim d^Z 1 / 2 u asJ = T x j„ zVj u aS)J} (48) 

on— shell 

where T x j Vs is calculated on-shell, from which the final cancellation of the gauge- dependence 
follows. 

If some of the flf do not vanish, the cancellations do not operate any longer and the 
amplitude turns out to be gauge parameter dependent f25(| . An explicit example has 
been considered in [fT5| |, for the W decay into quarks: if the CKM counterterm is gauge- 
dependent, the amplitude depends on the gauge parameters too. On the other hand, 
the above proof relies neither on a specific choice of renormalization of the unphysical 
parameters, nor on the regularization scheme adopted (provided the STI have been restored 
order by order). 



9 Summary 

We have introduced the Nielsen identities of the SM and used the problem of the definition 
of mass as a demonstrative example. In this context we have obtained some new results: 
we have proven to all orders in perturbation theory the gauge-parameter independence of 
the complex pole associated to any physical particle of the SM. We have considered the 
cases of the vector bosons, scalars and fermions in great generality, allowing for arbitrary 
mixing patterns. Particular attention has been paid to the case of the W boson, which is 
simpler because of the absence of mixing and has been chosen to illustrate some features 
common to all cases. Most of the proofs hold without modifications also in some extensions 
of the SM, like non-supersymmetric two-Higgs-doublet models. 

We have derived identities for the gauge-dependence of all the two-point functions 
of the SM, both for bosons and fermions, as well as for vertices and boxes involving 
external fermions. Using these expressions, we have shown the explicit mechanism of gauge 
cancellations which leads to gauge-independent four fermion amplitudes to all orders in 
the most general case of fermion and boson mixings and of CP violation. The formalism 
introduced in this paper, supplemented by the material given in App. A (the Lagrangian 
involving the BRST sources), should allow for a very simple derivation of the Nielsen 
identities for any proper Green function in the electroweak SM and in QCD. 

We have also extensively discussed the renormalization of the Nielsen identities with 
an arbitrary regularization, in the case the Nielsen identities (but not the STI) are broken 
by renormalization. In that case the identities are deformed by new terms, which we have 
identified in full generality and computed in a few cases of particular interest. We have 
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also derived new results concerning the infrared-finiteness of the W pole mass and the 
photon two-point function at q 2 = in the SM. For completeness, we report in App. C the 
expressions for the fermionic one- loop self-energies in a generic gauge. 

In conclusion, the formalism of the Nielsen identities can be useful in various appli- 
cations: (i) at the conceptual level, for the identification of gauge-independent quantities 
such as invariant charges and for the gauge-independent definition of renormalized pa- 
rameters [|15|]; (ii) at the practical level, because in higher orders calculations it is generally 
simpler to compute the gauge- dependence using the Nielsen identities, and because these 
identities allow for powerful checks. It deserves to be better known to theorists. 
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A. Nielsen identities for pedestrians 



The aim of this Appendix is to review very briefly the formalism of Slavnov- Taylor Iden- 
tities (STI) in the case of the Nielsen identities and to provide some material necessary 
for the explicit calculation of the Green functions involving the BRST sources. For a non- 



expert introduction to the STI for specific physical amplitudes, we refer to [19|. First, we 
recall that in our conventions the gauge-fixing term in the SM Lagrangian is given by 



GF 



2£ 4 



2£ 



c,w 



■A 2 ) 
fx 9 



(Al) 



We always set £ WjZ = £w'z = £w,z, i-e. we confine ourselves to the restricted 't Hooft 
gauge. Our starting point is the complete generating vertex functional T c , which generates 
the one-particle-irreducible Green functions. In order to simplify the structure of the 
STI, it is convenient to introduce for linear gauge-fixings a reduced generating functional 
r (sometimes indicated by T in the literature), which differs from T° by a local term, 
corresponding to the gauge-fixing part of the Lagrangian: 



r c - 



J d A x Cqf- 



(A2) 



In practice, the STI obtained from T coincide with the STI obtained from T c after imple- 
mentation of the ghost equation [p5 ]. Of course, one should keep in mind that the Green 
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functions involving unphysical fields generated by T coincide with the ones generated by 
r c only up to constant terms. For example, one has Ty) Wv = ^vPw v ~^P >1 P u /^w and 
r^| G _ = Tq+ g __ + £, W M^ at the tree level, while the difference at higher orders depends 
only on the renormalization of the W field and of the gauge parameters. As we have 
eliminated the classical gauge-fixing, it is clear that T w + G - ^ already at tree level. 

The invariance of the action under BRST transformations implies the STI for the 
functional T (see for ex. |^B| ) . 

s (r) = J d 4 x 

where and c° are the gauge boson for the U(l) abelian factor of the gauge group and 
the corresponding ghost, ip stands for any other quantum field in the SM Lagrangian 
(gauge fields, scalars, ghosts, and fermions) and 7^ is the BRST source associated to <p 
and is coupled to the non-linear BRST variation of (p in the classical action. In the case of 
a fermion // the spinorial source is denoted by r\i. We also introduce the Slavnov- Taylor 
operator Sr defined by 

By functional differentiation of Eq. ( |A3| ) with respect to some SM fields one gets the 
Slavnov- Taylor Identities (STI). Electric and ghost charge conservation, as well as Lorentz 
invariance, should be taken into account, according to the examples given in the text. 

In order to obtain the Nielsen identities for the gauge parameter dependence of irre- 
ducible Green functions, we have to consider the case of extended BRST symmetry |T|, 
which involves also the transformation of the gauge parameters; Eq. ( |A3|) takes then the 
form 

S(T) +]r Xl <%r = 0, (A5) 

i 

from which Eq. ([!]) follows after differentiating wrt \ and setting \ = 0. In the fermionic 
sector the expressions are slightly complicated by the anticommutation relations and the 
Nielsen identity becomes 

I 

where d x = djd\ and the arrows indicate the direction in which the functional derivative 
wrt the fermionic field acts (this is important for anticommuting fields). 

We have seen that both the Nielsen identities and the STI contain Green functions 
involving the BRST sources 7^ and rjf (for fermions) associated to the various fields of the 
SM. If we want to compute these Green functions at a given order in perturbation theory, 
we need to know how the sources are coupled to the fields. To this end, we give below 
the complete action involving the BRST sources, which can be useful as a reference and 



SB,, 



| ^ ST 



(A3) 



5j v Sep Sep #7^ 



(A4) 



r 5 5 d Y T d y r 5 6 r 



Sipi St]i 



5ipi 5t]j 



+ C0J <-> T)i) 



(A6) 
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to obtain the Feynman rules necessary for actual calculations involving 7^ and r/f. Apart 
from the well-known Feynman rules of the SM (see for instance the second paper in ) . 
nothing else is needed to evaluate the unconventional objects that appear in the identities. 
Using the convention = c w W^ + s w B^, where W^, are the third component of the 
triplet of SU(2) L and the Uy(l) gauge boson, respectively, we have 



C 



BRST 



u 

73 



[c w d^c z - s w df,c A - ig 
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c 



s w 



tec 



s w 



+ i a » {<V a - gJ abc Gy) - 7 C 3 {igc + c-} + 7cT {t| c d 
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d L + 9s—d L c a 



+ 
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^c z 
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-e Qu(c A + 



^c z 



w + 9 s —u u c a 



eQ d (c A + ) d K + 9 s-d H c a } + h.c. 



A a 



(A7) 



where A a are the Gell-Mann matrices, R and L indicate the right and left-handed compo- 
nents of the fermion fields, and s w = sxn9 w , c w = cos 6^. The hermitian conjugate for 
the fermionic part is added at the end. The ghost charge of the various sources, which is 
important in writing the STI, can be inferred by Eq. ( |A7| ), assigning a number +1 to the 
ghosts and requiring C to be ghost charge neutral. No BRST source needs to be introduced 
for the abelian vector field and for its ghost. 7^ is the source of the BRST transformation 
of the third component of the gauge boson triplet. 

The last ingredient for the calculation of the Green functions involving the source \i 
characteristic of the Nielsen identities, are the couplings of \ with the other fields. There 
is a source Xi associated to any gauge parameter £3 |. The relevant Lagrangian takes the 

5 Having set the two gauge parameters £.[ equal to each other, we can work with only one source Xi- 
This differs slightly from the procedure adopted in jl5|| , where two distinct sources X^ 1 ' 2 "* were kept. 
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form: 



2£ ; 

Xw 

2£w 



^ c a duG^ 



C 



(A8) 



B. Nielsen identities and regularization 



In this appendix we clarify the meaning of Eq. (|6]) and show how its structure is preserved 
if the STI are enforced at each perturbative order by means of appropriate non-invariant 
counterterms. 

Let us consider a non-invariant regularization, such as dimensional regularization in 
the implementation of Ref. [fT^l , and proceed to impose the renormalization conditions 
according to the procedure outlined in Sec. |[ At a given order n in perturbation theory 
the STI are violated. We now assume that at order n — 1 the STI have been restored by 
the introduction of appropriate non-invariant counterterms. Following the discussion of 
Sec. 0, the Nielsen identity corresponding to the extended BRST symmetry can be written, 
at order n — 1, in the following form (here we consider explicitly different gauge- fixing 
parameters) 



s(r) + £ x (%r + A Xi -r) 



(n-l) 



(Bl) 



where 



£ 

m=l 



.(n-l) 



Pij °t 



n-l 



E^(m) _ p(n-m-l) _ 
Xi 

m=l 



(B2) 



x 



5H{x) 



p(n— m— 1) 



The matrices ^ m \^^ m \8f' , and pff are straightforward extensions of the pa- 
rameters introduced in Eq. (|6|). Following the general theorem of renormalization theory 
known as Quantum Action Principle (QAP) f43fl , the terms breaking the Nielsen identity 
at order n are local polynomial of the fields and we have 



s(r) + 5>, (%r + A Xi -r) 



(n) 



AS B) + £XiAjW+0(xitt) 



(B3) 



where the new terms Aq + J2i Xi^xi are ^ oca ^ operators. We do not consider here the 
terms 0(xjXj) as they will not enter our forthcoming discussion. Now we can use the 
nilpotency of the operator Sr + J2i Xi% to establish the following consistency conditions 
for the breaking terms of Eq. (B3): 



<Sr Ao 



(n) 



0. 



S 



r 



A (n) 

A.2 



0. 



(B4) 
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In the absence of anomalies the first equation can be integrated obtaining the general 



solution [|ll],|i6 



A 



(n) 
CT 



(B5) 



where r^j, are local non-invariant counterterms. These counterterms are needed to restore 
the symmetries (in our case the STI) to the order n and are computed by standard tech- 

The removal of the breaking terms Aq™^ by 



niques of algebraic renormalization ||18|,|19| 
means of the counterterms I^y is essential in order to recover the unitarity of the theory 
and the physical interpretation of the S-matrix amplitudes. For what concerns the other 
breaking terms, namely , they do not play the essential role of the previous ones, but 

contain the information on the gauge dependence of T^. 

The new functional given by r^ n ^ + satisfies the STI identity at order n. On this 
basis we can study the gauge parameter dependence of the Green functions according to 
the Nielsen identities. Combining the second of Eqs. ( [B4]) with Eq. (|B5| ) we obtain 



0. 



(B6) 



where we have also used [d^, Sr ) = 0. Finally, the last equation can be solved using the 
cohomological methods outlined in Sec. 2, 



p(«) , A (n) 
CT ' 



X; 



X (n) +S To Y (n) . 



(B7) 



As discussed in the text, the terms in belong to the cohomology and represent the 
gauge parameter dependence of the physical parameters. On the other hand, the terms 
in Sr Y^ are cohomo logically trivial and contribute only to the unphysical parameters 
such as the renormalization of the fields, of the gauge fixing parameters etc. Therefore the 
insertion of the non-invariant counterterms at order n does not affect /3f'^ n \ 7^ n - ) , Sf 
at the same order and does not spoil the simple physical interpretation we have given them 
in the text. 

In summary, we have explicitly seen how the structure of Eq. (^]) is preserved at all or- 
ders. When the renormalization conditions are chosen according to the scheme presented 
in Sec. |2| and all the steps are properly performed, the result of the whole renormalization 
program are Green functions which at each order n are finite, satisfy the symmetry prop- 
erties of the model and provide S-matrix elements which are bound to be gauge-parameter 
independent. 



C. Gauge dependence of the fermionic self-energies 

In this appendix we present the explicit gauge-parameter dependence of the one-loop 
fermionic unrenormalized self-energies in the SM. We consider the most general case of 

6 In the framework of Ref. [Q the situation is more complicate as the operator Sr„ does not commute 
explicitly with the derivative wrt the gauge parameters. 
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mixing and define the fermionic self-energy Ejj as +i times the standard Feynman am- 
plitude for the transition j —*■ i and extract a factor g 2 . The expressions in the 't Hooft- 



Feynman gauge (£, = 1) can be found, for example, in Ref. [44]. At the one-loop level, 
instead of Eq. (|32[) , we can use the decomposition 



The individual components of the self-energies are then given in an arbitrary gauge by 
(similar formulae are also in ||45|| ) 



yS 
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5 



+ (G-l) 



„2 



n b zi +\e i r i i z M z + 



(CI) 



+ (£ 7 - 1) 5ij ~y Ql \p 2 {l - Xi) 2 c yi - (1 - Xi)a y - (1 + Xi)h 

+ (£ w - 1) ^2 ~r P 2 ( l ~ ?> x k)cwk - b wk - iwM w c wk - a, 
k 4 



2c 2 



- [£ 2 (l- Xt ) + - } n z 



X j 



(1 + X i)~J 



b Z i + iz M z c zi 



(C2) 



S 2 . r 

+ (£ 7 - 1) 5jj -y Q 2 [p 2 (l - Xj) 2 c 7i - (1 - Xj)a 7 - (1 + Xi)b ¥ 

Oiw - b W k + (ml + p 2 - £ w M w )c wk 



(U-i)5>5? 



4 p 2 



+ (6-1) 



2c?„ 



r 2 (l-x,) 2 



(1 + Xi 



r1{l + Xi) - 



1 i 



(b zi + £ z M 2 c z 



(C3) 



where we have used the following notation for the n-dimensional integrals (i, j = 7, Z°, W, f) 

d n k 1 



ft; 



b; ; 



A—n 



im 2 } 



4— n 



?,// 



4-n 



(2tt)" [fc 2 - m?][fc 2 - £. 

cTA; 1 
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We have also used X; 



m: 



''/p 2 , while 



If — Qi s w an d r i — —Qi s w are the l e ft an( i 



right-handed couplings of the fermion flavor i and Qi and If = ±5 its electric charge 
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and isospin. In the case of quarks, the mixing matrix factor A^ equals VikVA., where V 
is the CKM matrix, if i,j (k) are up (down) quarks and A^ = V^Vkj if the opposite is 
true. For leptons with massless neutrinos A^ = Sy 8^ Vi or 5^5^, i.e. there is no mixing. 
The gluon exchange diagrams can be obtained from the photonic ones setting Qi = 1 and 
multiplying by the color factor Cf- Notice that a 7 and c 7 j are infrared divergent and 
an infrared regulator (like a photon mass) should be introduced. Of course, the infrared 
divergences cancel out in Eqs.( P2| - |C3| ). It is straightforward to verify [(46| that in the 
diagonal case the mass counterterm, Smi/rrii = Sf^m^) + |S^(m^) + i£^(mf) +T i; where 
Tj is the tadpole contribution, is independent of the gauge parameters. From the off- 
diagonal parts of Eqs. (|Cl| - |C3l) it is easy to derive some of the results of Ref. ^5] on the 
gauge dependence of the CKM counterterm. 
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